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ABSTRACT 

An asymmetry between the ingress and egress times was observed in the near-UV light curve 
of the transit planet WASP- 12b. Such asymmetry led us to suggest that the early ingress in 
the UV light curve of WASP- 12b, compared to the optical observations, is caused by a shock 
around the planet, and that shocks should be a common feature in transiting systems. Here, 
we classify all the transiting systems known to date according to their potential for producing 
shocks that could cause observable light curve asymmetries. We found that 36/92 of known 
transiting systems would lie above a reasonable detection threshold and that the most promis- 
ing candidates to present shocks are: WASP-19b, WASP-4b, WASP-18b, CoRoT-7b, HAT-P- 
7b, CoRoT-lb, TrES-3, and WASP-5b. For prograde planets orbiting outside the co-rotation 
radius of fast rotating stars, the shock position, instead of being ahead of the planetary motion 
as in WASP- 12b, trails the planet. In this case, we predict that the light curve of the planet 
should present a late-egress asymmetry. We show that CoRoT-1 lb is a potential candidate to 
host such a behind shock and show a late egress. If observed, these asymmetries can provide 
constraints on planetary magnetic fields. For instance, for a planet that has a magnetic field 
intensity similar to Jupiter's field (~ 14 G) orbiting a star whose magnetic field is between 1 
and 100 G, the stand-off distance between the shock and the planet, which we take to be the 
size of the planet's magnetosphere, ranges from 1 to 40 planetary radii. 

Key words: planet-star interactions — planets and satellites: magnetic fields — stars: coronae 
— stars: winds, outflows 



1 INTRODUCTION 

Recently, Fossa tTet al.| j2010) observed that the near-UV transit 
light curve of the close-in giant planet WASP- 12b shows an early 
ingress as compared to its optical transit. Such observations were 
interpreted as due to the presence of asymmetries in the exosphere 
of the planet (e.g., |Lai et al.|20 10i. In particular, we showed that this 
asymmetry could be explained by the presence of a shock around 
the planet \ Vidotto, Jardine & Helling|2010| hereafter, VJH10). 

By analysing which characteristics of the ambient surrounding 
the planet could lead to an observable early UV ingress, we showed 
that shock formation is likely to be a common feature of transit- 
ing systems. Furthermore, once the stand-off distance of the shock 
(determined through the time difference of transit observations at 
different wavelengths) and the stellar magnetic field strength are 
known, we can derive the planetary magnetic field intensity. The 
planetary magnetic field is believed to be responsible for shielding 
the planet against the erosion of the planetary atmosphere by the 
host star's wind or the impact of energetic cosmic particles. Such 
effects could harm creation and development of life in the planet. 
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Furthermore, the presence of a planetary magnetic field may induce 
star-planet interactions, e.g., through reconnection between stellar 
and planetary magnetic field lines ( Stevens 2005 ). 

In this Letter, we extend the idea developed for WASP- 12b 
l |VJH10) to all known transiting exoplanets and investigate cases 
in which we expect to see light curve asymmetries. We discuss a 
selection of transiting planets for which we expect a shock to form 
(Section|2j. This shock can either cause an observable early ingress 
or late egress during transition (Section[3ji, depending on the char- 
acteristics of the ambient surrounding the planet, as well as on the 
planetary orbital properties. Conclusions are presented in Section]?] 



2 WHICH PLANETS MIGHT DEVELOP SHOCKS? 

As in |VJHl6"l we discuss two possible scenarios for shock forma- 
tion: (1) the stellar magnetic field is strong enough to confine the 
coronal plasma to co-rotate out to the orbital radius of the planet; 
(2) the stellar magnetic field is not sufficiently strong, and the coro- 
nal material escapes in the form of a stellar wind. The data used 
here were a selection of transiting systems from the exoplanets 
encyclopedia (http://exoplanet.eu/catalog.phpl, except for the sky- 
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projected stellar rotation velocity vsin(i'), which was taken from 
|Schlau fman (20101. In this Letter, we disregard planets that have 
either high eccentricity (> 0.3) or far orbits (> 0. 1 AU). 



2.1 Co-rotating hydrostatic corona 

For a shock to develop, the relative velocity between the planet 
and the surrounding stellar corona must be greater than the lo- 
cal sound speed ([VJH10]l: \u 9tCOt - Ug\ > c s . Here, u %cot = 
[v sin(i)/ sin(i)](R olb /R,) is the azimuthal velocity of the co-rotating 
corona at the orbital radius of the planet R or \,, c s = (k B T/m) i/2 
the sound speed, m = nm p the mean particle mass (we adopt 
fi = 0.66), m p the proton mass, k B the Boltzmann constant, and 
uk = (GM./Zfub) 1 ' 2 is the Keplerian velocity of the planet that 
orbits a star with mass M, and radius R,. Since the sound speed 
depends on the coronal temperature T, we can write this as the 
maximum temperature for shock formation 



(« v ,coi - uk) m/k B . 



(1) 



If the temperature of the local medium exceeds this value, then no 
shock will form. Table [JJ shows the maximum temperature for all 
the planets in our sample (an extended version of Table [JJ is avail- 
able online). We note that, for stars where we do not know the stel- 
lar rotation rate, we take the limiting case of zero stellar rotation 
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We note that the maximum temperature for shock formation for the 
majority of the planets is (1 - 4) x 10 6 K. Although stellar flares 
can raise the coronal plasma temperature above these values, it is 
unlikely that a corona would be hotter than this. 

For a shock to be detected, it must compress the local plasma 
to a density high enough to cause an observable level of optical 
depth. For a hydrostatic, isothermal corona, the local density is 
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where no is the density at the base of the corona. If we now con- 
sider T = T max , we can obtain a minimum density (column 12 of 
Table [TJ. Figure [JJ shows this critical density as a function of R olb 
for a range of planets, assuming that all stars have a base coro 
nal density equal to that of the Sun (n 0o = 10 8 cm 4 , Withbroe 
1988). For comparison, we show in Figure[JJthe minimum density 
(n b s = 1.5 x 10 6 cirT 3 ) derived for WASP-12b (upper dashed line, 
|VJH10) . The main uncertainties in the determination of the local 
density surrounding the planet come from the plasma abundances, 
cross-sections and geometric depths through the shock. The geo- 
metric effects alone allow for one order of magnitude uncertainty. 
We allow one further order of magnitude for the other effects. To 
conservatively account for these uncertainties in the density, we 
present in Figure[JJa lower limit for the density (lower dashed line) 
that could still provide detection (;j min ^ 10 4 cm 4 ). We note that 
36 out of 92 planets lie above a reasonable detection threshold, sug- 
gesting that a detectable shock might indeed be a common feature 
surrounding transiting planets. 

Some planets orbit extremely close to the co-rotation radius of 
the star. At this radius, u K ss u vfioc and Am ~ 0. Therefore, in these 
cases, namely, XO-3b, CoRoT-3b, CoRoT-4b, and CoRoT-6b, no 
shock will be developed. However, if these planets are embedded 
in the wind of their host stars, which does not co-rotate with the 
star, a shock may be formed. 
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Figure 1. Predicted coronal density at the orbit of each planet, assuming 
that the stellar magnetic field is strong enough to confine the coronal gas 
out to the orbit of the planet. We assume all the stars have the same solar 
coronal base density. The upper dashed line represents the minimum density 
value derived for WASP-12b (n obs = 1.5 X 10 6 cm 4 , |VJH10} . The lower 
dashed line represents a lower limit for the density given the uncertainties 
in abundances, cross-sections and geometric depths through the shock. Ap- 
proximately, 40% of the planets in our sample lie above the detection limit. 



Assuming that no = rco, Q probably represents effectively a 
lower limit on the base density because the Sun is a relatively in- 
active, slowly rotating star. For more rapidly rotating, active stars, 
the coronal density is likely to be higher. To estimate this effect, we 
scale n to the stellar angular velocity CI,, such that 



no = no, 



a 
a" 



v sin(i) 
2 km s~ 



(4) 



As Figure|2]shows, this has a modest effect on the overall range of 
predicted densities, but changes the order of planets with the high- 
est predicted coronal densities. Thus the three planets with the high- 
est predicted scaled densities are, in order: WASP- 18b, WASP- 19b, 
and WASP-12b. Figure[2]only contains the top 10 planets that may 
be observable with today's instrumentation, and the special case 
of CoRoT-llb (Table |2|. We note that, although the density pre- 
dicted for CoRoT-llb (right-most point in Figure [2]l is just above 
the threshold limit (lower dashed line), the assumption that the host 
star coronal base density should scale with CI, leads to a signifi- 
cant increase on the density of the ambient medium surrounding 
the planet, increasing the potential to detect light curve asymme- 
tries. We note that, in order for a transit to be clearly detected, 
the depth of the transit should be larger than the instrument's un- 
certainty. So, for instance, observations of the near-UV transit of 
CoRoT-7b, which possesses a relatively shallow transit, could be 
compromised if done with the Cosmic Origins Spectrograph on the 
HST (used in the case of WASP-12b). Observations of CoRoT-lb 
may also not result in a positive detection with such instrument, as 
its host star is about two magnitudes fainter (V band) than WASP- 
12. 

The stand-off distance of the shock from the planetary surface 
is also a factor that contributes to the observability of this effect. If 
we consider that the extent of the planetary magnetosphere r M can 
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Table 1. Transiting planetary systems known as of September/2010. Planets are ordered in increasing coronal densities of their host stars, taken as a proxy 
for the detection of light curves asymmetries. The columns are: (1) the planet name, (2) mass, (3) radius, (4) orbital period, and (5) semi-major axis, (6) the 
distance to the system, (7) the host star spectral type, (8) mass, and (9) radius, (10) the sky-projected stellar rotation velocity, (11) the maximum temperature 
required for shock formation; (12) the local density around the planet for the confined corona, and (13) considering the coronal density scales with fi,, (14) 
the size of the planet magnetosphere for a planet with B p = 14 G and a star with B, = 1 G, (15) the same but for B„ = 100 G, (16) the minimum planetary 
magnetic field relative to the stellar one (/ = (B p / B t ) m ] n ) that is required to sustain a magnetosphere. An extended version of this table is available online. 
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Figure 2. Influence of different assumptions on the predicted density at the 
orbit of each planet in Table|2] for the co-rotating corona with base density 
that is unsealed (squares) and scaled with fl, (crosses) densities, and for the 
wind case (triangles). Dashed lines are as in Figure[T] 



be taken as approximately the observationally derived stand-off dis- 
tance of the shock, we can calculate the magnetic field intensity of 
the planet ( VJH10). In the case where the shock location is deter- 
mined simply by static pressure balance, we have 



mn c Aw B c (R olb ) 

+ : + Pc 



B„(r M ) 2 



(5) 



2 4n r " An 

where n c , p c and B c (R olb ) are the local coronal density, thermal 
pressure, and magnetic field intensity, and p p and B p (r M ) are 
the planet thermal pressure and magnetic field intensity at r M . 
In |VJH10l we showed for WASP- 12b that the dominant terms in 



Equation |5| are the magnetic terms. We assume the same is true in 
our next calculations. Therefore, considering that both the star and 
the planet magnetic field geometry can be approximated as that of a 
dipolar field, we have B,(R OTb /R,)~ 3 ^ B p (r M /R p y 3 , where B, and 
B p are the magnetic field strength at the surfaces of the star and the 
planet, respectively. Hence 
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Therefore, for a given ratio B p jB„ r M scales linearly with R olb . 
Figure [3] shows the effect on r u due to a change in B,. The planets 
are all assumed to have B p = 14 G, similar to that of Jupiter. We 
obtain values for r M ranging from 1 to 40 R p for B, = 1 to 100 G. 
These values are comparable to those for the solar system planets. 
If we take the stand-off distance between the shock and the planet 
to be r M , planets with large r M present earlier ingresses (or later 
egresses) than planets with small r M - 

Not all planets may be capable of generating a magnetic field 
of sufficient strength to sustain a magnetosphere. Table[T|shows the 
values of the minimum ratio / = (B p /B,) min required to support a 
magnetosphere. Note that a planetary magnetic field of about a few 
percent of that of the host star is enough to support the existence 
of a planetary magnetosphere. If the kinetic and thermal terms in 
Equation [5] are important, the minimum B p would even be smaller. 



2.2 Stellar Wind 

If the coronal material is not magnetically confined, it escapes in 
the form of a wind. By assuming that the host star has an isother- 
mal, thermally-driven wind (Parker 1958), we calculate the wind 
velocity and density profiles for a few selected planetary systems. 
Systems that were top-ranked on Table[T]but with distances > 1 kpc 
where excluded from our calculations, because they are too far 
away to be observable with the instrumentation available today. 
One special case (CoRoT-1 lb) will be discussed in Section[3] 
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Figure 3. Magnetospheric size in units of the planetary radius for a stellar 
magnetic field strength of 1 G (triangles) and 100 G (diamonds). All planets 
are assumed to have a field of 14 G. Planets that are surrounded by coronal 
material denser than 10 4 cm~ 3 (detection limit) are represented by filled 
symbols. If we take the stand-off distance between the shock and the planet 
to be the size of the planet's magnetosphere, planets with large tm present 
earlier ingresses (or later egresses) than planets with small ru- 



The wind radial velocity u r is derived from the integration of 
the differential equation 



pu r 



du r 



dp GM, 
dr ^ r 2 



(7) 



where p = nm is the wind mass density, r is the radial distance from 
the star, p = nk B T mxi is the thermal pressure. We assume that the 
temperature of the wind is the maximum temperature for shock for- 
mation given in Table p[f| Table[2]presents the wind characteristics 
at R olb for the selected planetary systems. For the top-ranked plan- 
ets, the u r ranges from 240 - 360 km s -1 (column 8). The derived 
local densities (column 11) are obtained considering that the wind 
base density is the solar value. They are also shown in Figure [2] 
(filled triangles) and are 2.5 - 3.5 times smaller than the values de- 
rived from the confined corona scenario (shown in column 12 for 
comparison). The derived mass-loss rates of the host stars (column 
14) are at most two orders of magnitude larger than the solar value. 

Table [2] also presents the azimuthal velocities obtained for 
the wind and corona scenarios (columns 9 and 10, respectively). 
Along with the Keplerian velocity u K of the planet, they provide 
the condition for shock formation: \u K — M ViCor | > c s for the corona 
scenario, and [\u K - ^ >w j n d| 2 + k^] 1/2 > c s for the wind scenario. 
Note that if the corona is confined by the stellar magnetic field, 
and thus co-rotates with the star, w v>cor = [vsin(i)/ sm(i)](R olb /R,). 
For the wind case, the azimuthal velocity is determined by con- 
servation of angular momentum of the particles leaving the star 

^^.wind 

= [vsin(0/ sin(0](7?*/i?orb)- For the calculations in this pa- 
per, we adopt sin(;) = 1. If, however, sin(;) t 1, then the actual 



Note that in the wind scenario, the maximum temperature is actually 
given by the condition that An 2 + u 2 ■ 

7"max,wuid = [(" v ,wind - U K ) 2 + uj]m/k B 



The angle that the shock normal makes to the relative az- 
imuthal velocity of the planet is 



arctan 
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("v.wind -Uk) + u- > c-. Therefore, 



\ \UK — «<t,windl 

Such values are also presented in Table [2] (column 13), where we 
note that 9 ranges between 48 - 58°. As shown in VJH10| 6 in- 
creases with the wind temperature, due to the increase caused in u r 
(the thermally-driven wind is more accelerated for large T). 



3 LATE EGRESS OR EARLY INGRESS? 

For the ahead-shock case, which occurs in the confined corona sce- 
nario, 9 = 0°, and the shock is always formed perpendicular to the 
line that connects the centres of the planet and the star. Another lim- 
iting case, which was only briefly discussed in VJH10, is the case 
when the shock trails the planet {8 = 180°). In the "behind-shock" 
case, the planet orbits the star beyond the Keplerian co-rotation ra- 
dius, so that the coronal plasma lags behind the planetary motion. 
In other words, for a circular orbit, u K - u ffioI < or, similarly, 
P, < .P or b, where P, = 2^/?,/[vsin(0/ sin(i)] is the stellar rota- 
tional period. In order to develop a behind-shock, the planet must 
be in a prograde orbit. 

It is interesting to note that, in the behind-shock limit, the 
asymmetry in the time of ingress/egress of the light curve will be 
different from the one observed in WASP- 12b, where an ahead- 
shock is expected to be present ( VJHIOl. Because the compressed 
material trails the planet, in the behind-shock case, the light curve 
asymmetry occurs at the time of the egress: the planet may have 
completed the transit through the stellar disk, but the compressed 
material will still be seen in "transit" for an extra amount of time. 
During this time, if the material is able to absorb stellar radiation 
in the Mg II lines, for instance, the transit will appear to us as a 
late egress instead of an early ingress, as compared with WASP- 
12b. Planets with a certain mass escaping rate can also develop 
cometary tails due to evaporating processes caused by the ionis- 
ing radiation emitted by the star (Ehrenreich et al. 2008). Because 
a cometary tail around the planet would also produce a late egress, 
the behind-shock scenario can hence mimic a cometary tail sce- 
nario. If CoRoT-7b indeed hosts a cometary tail, as suggested by 
|Mura et aL] l |2010| l, then its light curve asymmetry, if observed, will 
show a late egress. However, we predict the existence of an early 
ingress caused by the presence of an ahead-shock. 

The condition for a behind shock (P , < P mD ) was met in 5 of 
the 71 cases where measurements of v sin(i) exist, namely: WASP- 
33b, CoRoT-1 lb, WASP-7b, CoRoT-4b, and CoRoT-6b. WASP-7b, 
CoRoT-4b, and CoRoT-6b present 10 R, < R orb < 18 R,, which 
may be too distant to consider that the corona co-rotates with the 
star. In these three cases, the impacting wind is expected to signifi- 
cantly affect the geometry of the shock, which should be formed at 
an intermediate angle 9. WASP-33b is an interesting case to study 
in future works, because its host star is in another mass regime. 
Therefore, it may have a different internal structure, it should host 
a different magnetic field structure, which may not confine the stel- 
lar material up to the orbit of the planet. As a consequence, the best 
candidate to host a behind-shock among our sample of planets and, 
therefore, present a late-egress in the UV light curve asymmetry, is 
CoRoT-llb. 

CoRoT-1 lb is an example of a prograde giant gas planet or- 
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Table 2. Wind calculation for selected planets. The columns are: (1) the planet name, (2) the stellar mass, (3) and radius, (4) the planet semi-major axis, (5) 
the maximum temperature necessary for shock formation, (6) the sky-projected stellar rotation velocity, (7) the Keplerian velocity of the planet, (8) the radial 
velocity of the wind at R or b, (9) and (10) the azimuthal velocities of the wind and of the corona at the planetary orbital distance, (1 1) and (12) the density at 
the planetary orbital distance for the isothermal wind model and for the coronal model, (13) the angle 6 of the shock normal, and (14) the wind mass-loss rate. 
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biting a rapidly rotating F6-star (vsin(i) = 40 km s , Gandolfi 



|et al.||20"To| >. If the magnetic field of CoRoT-11 is strong enough 
to maintain the coronal plasma in co-rotation up to the orbital ra- 
dius of the planet (R OI b = 6.84 R,), then the azimuthal velocity of 
the corona at the planet orbit is w ViCor - 274 km s . With an or- 
bital period of P olb = 2.99 d, CoRoT-1 lb has a Keplerian velocity 
of u K = 161 km s -1 , which leads to a relative azimuthal velocity 
Am = \tig — w™ lCOT | = 113 km s -1 . In order for a shock to form behind 
the planet, T < T max =s 10 6 K (Equation (]}). 

If the coronal plasma is not co-rotating with the star at the or- 
bit of CoRoT-llb, the host star's wind influences shock formation 
around the planet. By assuming an isothermal wind at 10 6 K, the 
shock normal is formed at 9 = 36° (Table [2]l. However, the low 
ambient density around the planet (~ 25 times smaller in the wind 
model compared to the confined corona model), may present some 
problems for the detection of the absorption during the UV transit. 
Such a small density is reflected in the sub-solar wind mass-loss 
rate. It is interesting to note that, although the planet is located in 
the sub-sonic region of its host star wind, a shock is formed because 
of the planet's high orbital velocity. 

Therefore, through observations of the planet light curve 
asymmetry (late egress/early ingress), one can constrain whether 
CoRoT-1 1 has a confined corona up to R mb (6 = 180°, late egress), 
or its wind affects shock formation (8 =s 36°, early ingress), provid- 
ing that the coronal temperature allows shock formation. 



4 CONCLUSION 

In this Letter, we determined which planets are prone to develop a 
surrounding shock that could be observed in asymmetries of plan- 
etary light curves. We discussed two possible scenarios for the for- 
mation of a shock: (1) when the stellar magnetic field is strong 
enough to confine the hot coronal plasma out to the planetary or- 
bit and (2) when the stellar magnetic field is unable to confine the 
plasma, which escapes in a wind. For both cases, we predicted the 
characteristics of the ambient medium that surrounds the planet, 
and discussed whether such characteristics presents favourable con- 
ditions for the presence of a shock. The planets that were top ranked 
according to scenario (1) are: WASP-12b, OGLE-TR-56b, WASP- 



19b, SWEEPS- 11, WASP-4b, WASP-18b, CoRoT-7b, CoRoT-14b, 
HAT-P-7b, OGLE-TR-132b, CoRoT-lb, TrES-3, and WASP-5b. 

With the observational determination of the stand-off distance 
between the shock and the planet (which we consider to be the 
size of the planet's magnetosphere r M ), one can infer the planetary 
magnetic field intensity, provided the stellar magnetic field inten- 
sity is known. By assuming that the planet hosts a magnetic field of 
the same intensity as Jupiter's and the star's magnetic field is in the 
range between 1 and 100 G, we showed that r M ^ 1 - 40 R p . In gen- 
eral, for the stars in our sample, we show that, to sustain a planetary 
magnetosphere, a minimum planetary field intensity is required to 
be only a few percent of the stellar magnetic field. 

For the cases where the star is sufficiently magnetised to con- 
fine the coronal material, the presence of a shock ahead of the plan- 
etary motion can cause an early ingress in the light curve of the 
planet, while a shock that trails the planet (behind shock) can re- 
sult in a late egress. We showed that prograde planets orbiting fast- 
rotating stars can develop a behind-shock. In particular, CoRoT-1 lb 
is a promising candidate to host a behind shock. 
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6 Vidotto, Jardine, & Helling 



Table 1. Online material. Transiting planetary systems known as of September/2010. Planets are ordered in increasing coronal densities of their host stars, 
taken as a proxy for the detection of light curves asymmetries. The columns are: (1) the planet name, (2) mass, (3) radius, (4) orbital period, and (5) semi-major 
axis, (6) the distance to the system, (7) the host star spectral type, (8) mass, and (9) radius, (10) the sky-projected stellar rotation velocity, (1 1) the maximum 
temperature required for shock formation; (12) the local density around the planet for the confined corona, and (13) considering the coronal density scales 
with fi„, (14) the size of the planet magnetosphere for a planet with B p = 14 G and a star with S, = 1 G, (15) the same but for B, = 100 G, (16) the minimum 
planetary magnetic field relative to the stellar one (/ = (B p /B t ) mm ) that is required to sustain a magnetosphere. 
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F 


1.38 


1.81 


8.5 


0.87 


3.19 


3.56 


14.5 


3.1 


0.45 


Lupus-TR-3b 


0.81 


0.89 


3.91 


0.046 




K1V 


0.87 


0.82 




1.34 


3.17 




29.2 


6.3 


0.06 


HAT-P-19b 


0.29 


1.13 


4.01 


0.047 


215 


K 


0.84 


0.82 




1.29 


3.15 




29.3 


6.3 


0.05 


WASP-28b 


0.91 


1.12 


3.41 


0.046 


334 


F8-G0 


1.08 


1.05 


2.2 


1.25 


3.14 


3.16 


22.4 


4.8 


0.12 


HD 189733b 


1.15 


1.15 


2.22 


0.031 


19.3 


K1-K2 


0.80 


0.79 


3.32 


1.20 


3.04 


3.36 


20.6 


4.4 


0.16 


WASP- 11 
HAT-P-lOb 


0.46 


1.05 


3.72 


0.044 


125 


K3V 


0.82 


0.81 


0.5 


1.22 


2.95 


2.44 


28.0 


6.0 


0.06 


HAT-P-12b 


0.21 


0.96 


3.21 


0.038 


142.5 




0.73 


0.70 


0.5 


1.24 


2.88 


2.43 


28.3 


6.1 


0.06 


WASP-21b 


0.30 


1.07 


4.32 


0.052 


230 


G3V 


1.01 


1.06 


1.5 


1.07 


2.67 


2.52 


25.3 


5.5 


0.09 
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J 
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^j^rn/ s f 


(MK) 
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(l) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


CO 


(12) 


(13) 


(14) 


(15) 


(16) 


XO-lb 


0.90 


1.18 


3.94 


0.049 


200 


G1V 


1.00 


0.93 


1.11 


1.21 


2.59 


2.37 


27.1 


5.9 


0.07 


WASP-22b 


0.56 


1.12 


3.53 


0.05 


300 




1.10 


1.13 


3.5 


1.04 


2.47 


2.66 


21.4 


4.6 


0.14 


WASP-2b 


0.85 


1.04 


2.15 


0.031 


144 


K1V 


0.84 


0.83 


5 


1.04 


2.41 


2.88 


19.4 


4.2 


0.19 


WASP- 16b 


0.86 


1.01 


3.12 


0.042 




G3V 


1.02 


0.95 


3 


1.12 


2.30 


2.50 


23.0 


5.0 


0.11 


WASP-29b 


0.24 


0.79 


3.92 


0.046 


80 


K4V 


0.82 


0.85 


1.5 


0.96 


1.84 


1.79 


27.9 


6.0 


0.06 


HAT-P-18b 


0.20 


1.00 


5.51 


0.056 


166 


K 


0.77 


0.75 




0.98 


1.49 




38.5 


8.3 


0.02 


OGLE-TR-lOb 


0.68 


1.72 


3.10 


0.042 


1500 


GorK 


1.18 


1.16 


7 


0.88 


1.41 


1.89 


18.5 


4.0 


0.22 


HD 209458b 


0.64 


1.38 


3.52 


0.047 


47 


GOV 


1.00 


1.15 


4.7 


0.73 


0.96 


1.28 


21.4 


4.6 


0.14 


WASP-25b 


0.58 


1.26 


3.76 


0.047 


169 


G4 


1.00 


0.95 


3 


0.88 


0.85 


1.04 


25.7 


5.6 


0.08 


CoRoT-8b 


0.22 


0.57 


6.21 


0.06 


380 


K1V 


0.88 


0.77 




1 


0.83 




42.2 


9.1 


0.02 


HAT-P-15b 


1.95 


1.07 


10.86 


0.096 


190 


G5 


1.01 


1.08 




0.75 


0.13 




46 A 


10.0 


0.01 


CoRoT-2b 


3.31 


1.47 


1.74 


0.03 


300 


G7V 


0.97 


0.90 


11.85 


0.74 


-0.05 


0.76 


16.1 


3.5 


0.33 


HAT-P-lb 


0.52 


1.22 


4.47 


0.055 


139 


GOV 


1.13 


1.12 


3.75 


0.73 


-0.40 


-0.17 


25.6 


5.5 


0.08 


HAT-P-8b 


1.52 


1.50 


3.08 


0.049 


230 




1.28 


1.58 


11.5 


0.48 


-1.42 


-0.85 


15.9 


3.4 


0.34 


HAT-P-llb 


0.08 


0.45 


4.89 


0.05 


38 


K4 


0.81 


0.75 


1.5 


0.71 


-1.46 


-1.46 


36.5 


7.9 


0.03 


HAT-P-6b 


1.06 


1.33 


3.85 


0.052 


200 


F 


1.29 


1.46 


8.7 


0.53 


-1.51 


-1.04 


18.5 


4.0 


0.22 


Kepler-8b 


0.6 


1.42 


3.52 


0.05 


1330 




1.21 


1.49 


10.5 


0.47 


-1.75 


-1.20 


16.8 


3.6 


0.29 


XO-4b 


1.72 


1.34 


4.13 


0.06 


293 


F5V 


1.32 


1.55 


8.8 


0.49 


-1.94 


-1.49 


18.5 


4.0 


0.22 


GJ 436b 


0.07 


0.37 


2.64 


0.029 


10.2 


M2.5 


0.45 


0.46 


2.4 


0.59 


-2.05 


-1.64 


32.1 


6.9 


0.04 


GJ 1214b 


0.02 


0.24 


1.58 


0.014 


13 




0.16 


0.21 


2 


0.41 


-3.18 


-2.51 


34.2 


7.4 


0.03 


WASP- 17b 


0.49 


1.74 


3.74 


0.051 




F6 


1.20 


1.38 


9 


0.43 


-3.48 


-2.97 


19.1 


4.1 


0.20 


HAT-P-14b 


2.2 


1.2 


4.63 


0.06 


205 


F 


1.39 


1.47 


8.4 


0.41 


-5.42 


-4.96 


20.9 


4.5 


0.15 


WASP- 1 Ob 


3.06 


1.08 


3.09 


0.037 


90 


K5 


0.71 


0.78 


6 


0.39 


-5.83 


-5.25 


24.4 


5.3 


0.09 


OGLE-TR-lllb 


0.54 


1.08 


4.01 


0.047 


1500 


GorK 


0.82 


0.83 


5 


0.33 


-10.13 


-9.65 


29.2 


6.3 


0.06 


WASP-7b 


0.96 


0.92 


4.95 


0.062 


140 


F5V 


1.28 


1.24 


17 


0.17 


-24.53 


-23.70 


25.8 


5.6 


0.08 


HAT-P-9b 


0.78 


1.40 


3.92 


0.053 


480 


F 


1.28 


1.32 


11.9 


0.16 


-27.30 


-26.64 


20.7 


4.5 


0.16 


HAT-P-23b 


2.09 


1.37 


1.21 


0.023 
























HAT-P-22b 


2.15 


1.08 


3.21 


0.041 
























WASP-6b 


0.50 


1.22 


3.36 


0.042 


307 


G8 






1.4 












0.00 


XO-3b 


11.79 


1.22 


3.19 


0.045 


260 


F5V 


1.21 


1.38 


18.54 


0.04 






17.0 


3.7 


0.28 


WASP-36b 


2.40 


1.40 


1.50 


























HAT-P-20b 


7.25 


0.87 


2.88 


0.036 
























HAT-P-21b 


4.06 


1.02 


4.12 


0.049 
























CoRoT-3b 


21.66 


1.01 


4.26 


0.057 


680 


F3V 


1.37 


1.56 


17 


0.01 






18.9 


4.1 


0.21 


CoRoT-4b 


0.72 


1.19 


9.20 


0.090 




F8V 


1.10 


1.15 


6.4 


0.00 






40.4 


8.7 


0.02 


CoRoT-6b 


2.96 


1.17 


8.89 


0.086 




F9V 


1.06 


1.03 


7.5 


0.07 






43.0 


9.3 


0.02 


OGLE2-TR-L9b 


4.34 


1.61 


2.49 




900 


F3 


1.52 


1.53 


39.33 














WASP- 13b 


0.46 


1.21 


4.35 


0.053 
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G1V 






4.9 














WASP-31b 


0.50 


1.60 


3.50 



























